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ABSTRACT 

We investigate the axisymmetric magnetosphere of an aligned rotating mag- 
netic dipole surrounded by an ideal force-free plasma. We concentrate on the 
magnetic field structure around the point of intersection of the separatrix be- 
tween the open and closed field-line regions and the equatorial plane. We first 
study the case where this intersection point is located at the Light Cylinder. 
We find that in this case the separatrix equilibrium condition implies that all the 
poloidal current must return to the pulsar in the open-field region, i.e., that there 
should be no finite current carried by the separatrix/equator current sheet. We 
then perform an asymptotic analysis of the pulsar equation near the intersection 
point and find a unique self-similar solution; however, a Light Surface inevitably 
emerges right outside the Light Cylinder. We then perform a similar analysis for 
the situation where the intersection point lies somewhere inside the Light Cylin- 
der, in which case a finite current flowing along the separatrix and the equator is 
allowed. We find a very simple behavior in this case, characterized by a 90-degree 
angle between the separatrix and the equator and by finite vertical field in the 
closed-field region. Finally, we discuss the implications of our results for global 
numerical studies of pulsar magnetospheres. 
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1. Introduction 

The structure of the pulsar magnetosphere has been an active area of research for more 
that 30 years, starting with the pioneering work by Goldreich & Julian (1969). Despite 
the fact that real astrophysical radio-pulsars are believed to be oblique rotating magnetic 
dipoles, much of the theoretical effort has been devoted to a significantly simpler case of an 
aligned rotating magnetic dipole (i.e., a "non-pulsing pulsar"), in which case the problem 
becomes stationary and axisymmetric. The situation becomes even more tractable if one 
also assumes the space around the pulsar to be filled with a plasma that, on one hand, is 
dense enough to shorten out the longitudinal electric fields, thus providing the basis for the 
ideal-magnetohydrodynamic approximation, and, on the other hand, is at the same time 
tenuous enough for all the non-electromagnetic forces to be negligible, thus enabling one to 
regard the plasma as force- free. 

All these simplifications have lead early on to the derivation, simultaneously by several re- 
searchers, of the main equation governing the structure of the magnetosphere, the so-called 
pulsar equation (Scharlemann & Wagoner 1973; Michel 1973b; Okamoto 1974), which is 
essentially a special-relativistic generalization of the well-known force-free Grad-Shafranov 
equation. This equation is a quasi-linear elliptic second-order partial differential equation 
with a regular singular surface, the so-called Light Cylinder. Despite the fact that this, in 
general non-linear, equation was first derived thirty years ago, most of the attempts to solve 
it, both analytical and numerical, have been limited, until very recently, to its linear spe- 
cial cases and to the region inside the Light Cylinder (e.g., Scharlemann & Wagoner 1973; 
Michel 1973a; Beskin et al. 1983; Beskin & Malyshkin 1998). These self-imposed restrictions 
can of course be attributed to the considerations of mathematical convenience and tractabil- 
ity; however, they are not very well physically motivated and thus may be too restrictive to 
be relevant to real astrophysical systems. 

At the same time, by the late 1990s, the available computing resources and numerical tech- 
niques have become sufficiently powerful for the most general nonlinear problem, involving 
the regions on both sides of the Light Cylinder, to become tractable. In their pioneering work, 
Contopoulos et al. (1999, hereafter CKF) have obtained the first (and apparently unique) 
numerical solution of the general problem, albeit with a rather poor spatial resolution. Their 
approach has later been used by Ogura & Kojima (2003) who have obtained essentially the 
same results, but with a higher numerical resolution. Both of these groups have focussed on 
the global structure of the solution; consequently, they have not paid enough attention, in 
our opinion, to some key but subtle issues regarding the separatrix between the open and 
closed field-line regions and, especially, the point of intersection of this separatrix with the 
equatorial plane. As a result, the separatrix has not, we believe, been treated correctly, and, 
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in particular, the separatrix equilibrium condition (Okamoto 1974; Lyubarskii 1990) has not 
been satisfied close to the intersection point. We thus suspect that their solution is probably 
not quite correct. 

In this paper we try to clarify how the magnetic field near the separatrix and especially 
around the separatrix-equator intersection point should be treated, as this question appears 
to be extremely important for setting up the correct boundary conditions for the global 
problem. 

In § 2 we review the basic equations describing the ideal-MHD force- free magnetosphere of an 
aligned rotating dipole; we pay special attention to the so-called Light Cylinder regularity 
condition and discuss how this condition could be used to find the correct form of the 
poloidal-current function and, simultaneously, to obtain the solution that passes smoothly 
across the Light Cylinder. Then, in § 3, we investigate some general properties of the 
magnetic field near the putative separatrix between the open and closed field regions; we use 
the separatrix equilibrium condition to show that, if the separatrix intersects the equator 
at the Light Cylinder (as has been assumed in the numerical simulations by CKF and by 
Ogura & Kojima 2003), then the poloidal current and hence the toroidal magnetic field have 
to vanish on the last open field line above the separatrix. We then use these findings in § 4 to 
construct a unique self-similar asymptotic solution of the pulsar equation in the vicinity of 
such an intersection point; in particular, we find all the power-law exponents describing the 
field near this point and also the angle between the separatrix and the equator. At the end of 
this section we show that a Light Surface (a surface where the electric field becomes equal to 
the magnetic field and therefore where the particle drift velocity reaches the speed of light) 
has to appear right outside the Light Cylinder (it starts at the intersection point and extends 
outward at a finite angle with respect to the equator). We thus conclude that the only way 
to get a force-free solution that would be valid at least some finite distance beyond the Light 
Cylinder, is to consider the case of a separatrix-equator intersection point lying some finite 
distance inside the Light Cylinder. We investigate the magnetic field structure around such 
a point in § 5. We find that in this case the poloidal-current function may stay finite, but 
its derivative with respect to the poloidal flux has to go to zero on the last open field line; 
we then find the asymptotic behavior of the magnetic field around the intersection point 
corresponding to this case. Finally, we present our conclusions in § 6, where we also discuss 
the implications of our results for the past and future numerical studies of the axisymmetric 
pulsar magnetospheres. 
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2. Basic Equations 

We consider the magnetosphere of an aligned rotating magnetic dipole under the assumptions 
of stationarity, axial symmetry with respect to the rotation axis Z, and reflection symmetry 
with respect to the equator Z = 0. We shall include the effects of special relativity but will 
ignore general- relativistic effects (i.e., work in Euclidean space). We shall assume that the 
dipole is surrounded by a very tenuous, but highly conducting plasma. The first of these 
assumptions enables us to neglect all non-electromagnetic forces (i.e., gravity, inertial, and 
pressure forces), and thus to conclude that the structure of the magnetosphere is governed 
by the relativistic force-free equation: 

pE + ^-^ = 0. (i) 

Here, the electric charge density p and the electric current density j are related to the electric 
and magnetic field through the steady-state Maxwell equations: 

p = -Lv-E, (2) 
j = ^VxB. (3) 

The second of our assumptions concerning the magnetospheric plasma, i.e., the assumption 
of infinite conductivity, enables us to use ideal magnetohydrodynamics (MHD): 

E+^ = 0. (4) 

c 

These equations (with the appropriate boundary conditions) should be sufficient for deter- 
mining the structure of the magnetosphere. 

As is well known, an axisymmetric magnetic field B can be described in terms of two func- 
tions, the poloidal magnetic flux (per one radian in the azimuthal direction) ^ and the 
poloidal current /, as 

B = WxV0 + /V0, (5) 

where is the azimuthal (or toroidal) angle. In cylindrical coordinates (R, 0, Z) the magnetic 
field components then are: 

= -R8Z> < 6 > 
= ^- (8) 
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Next, by applying Faraday's law dB/dt = — cV x E in a steady state, we see that V x E = 0, 
and hence E = — V$; together with the assumption of axial symmetry, this gives 

£0 = 0, (9) 

that is, the electric field is purely poloidal. Since, as it follows from equation (1) or (4), the 
electric field must be perpendicular to the magnetic field, we can write it as 

E = -B pol x^ = --Vf, (10) 

c c 

where B pol = x V0 is the poloidal magnetic field, and is the unit vector in the 0- 
direction. By substituting this relationship into V x E = 0, we find that f2, which has 
the meaning of the angular velocity of magnetic field lines, is constant along the field lines, 
f2 = This is the well-known Ferraro isorotation law. 

In the problem we are interested in, all the field lines are tied, at least at one end, to the 
surface of the pulsar, which is assumed to be rotating as a solid body with some uniform 
angular velocity We also assume that our ideal-MHD assumption is valid all the way 
up to the pulsar surface, i.e., that there is no appreciable gap between this surface and the 
ideal-MHD region. Therefore, from now on, we shall assume that all the field lines rotate 
with the same angular velocity 2 

n(W) = ft* = const . (11) 



From equation (10) we find 

|E| = |B pol |x, (12) 

where 

x = - = rTc (13) 

is the cylindrical radius normalized to the radius of the Light Cylinder (LC) R LC = c/Q. 
Thus, we see that |E| < |B po i| inside the LC and |E| > |B po i| outside the LC. 

Now let us turn to the force-free equation (1). The toroidal component of this equation, 
together with equation (9), gives: 

jpoi x B pol = =s> [V/ x V0] x [W x V0] = =s> / = , (14) 



2 The case where 0(4') ^ f2* was considered, e.g., by Beskin et al. (1983) and by Beskin & Malyshkin 
(1998). 
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and the poloidal component (perpendicular to B po i) of equation (1) can be written as 



d 2 ^ d 2 ^ 



1 + x 2 d^I 





1 


5k = 


xR 2 LC dz 




1 


B z = 


xR 2 LC dx 


Bq = 


m 


xR 2 LC 



where we have introduced z = Z/Rlc and where we also incorporated Rlc into /(^IO on the 
right-hand side (RHS) by defining 

/(*) = R L cI{^) • (16) 
Then the magnetic field components can be rewritten in terms of ^(x, z) as 

(17) 
(18) 
(19) 

Equation (15) is the famous pulsar equation (Scharlemann & Wagoner 1973; Michel 1973b; 
Okamoto 1974), also known in literature as the relativistic force-free Grad-Shafranov equa- 
tion (e.g., Beskin 1997). It is an elliptic second-order partial differential equation (PDE) for 
the flux function z); the left-hand side (LHS) of this equation is linear, whereas the 
RHS is, in general, nonlinear. 

One very important feature of equation (15) is that it has a regular singular surface at the 
Light Cylinder x — 1. Both indicial (or characteristic) exponents for this equation are equal 
to zero and so a general solution of equation (15) diverges logarithmically at this surface 
(e.g., Bender & Orszag 1978). On physical grounds, one imposes an additional condition 
that both the function ty(x, z) and its 1st and 2nd derivatives remain finite near the LC. For 
any given function it is indeed possible to find such non-divergent solutions separately 

on each side of the LC. Then, as can be seen from equation (15) itself, the contribution of the 
first term on the LHS [the term proportional to (1 — x 2 )] vanishes as x — > 1, and thus such 
non-divergent solutions on both sides of the LC have to satisfy the well-known LC regularity 
condition (e.g., Scharlemann & Wagoner 1973; Okamoto 1974; Beskin 1997): 

^(x = l,z) = lim. (20) 



This condition is very important and merits a few extra words of discussion. Our force-free 
equation with a constant field-line angular velocity Q is a 2nd-order equation with one, a 
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priori unknown, integral of motion, I (iff), and one singular surface, the LC. According to 
general theory (Beskin 1997), one then needs two boundary conditions to set up the problem 
properly. We believe that the correct approach here is to set the boundary conditions for 
the function iff at both the inner boundary (inside the LC: on the surface of the star, the 
symmetry axis, the separatrix between the open and closed field-line regions, and maybe the 
equator) and at whatever outer boundary (outside the LC: the equator and infinity or the 
Light Surface). At the same time the poloidal current function I (iff) should not be prescribed 
explicitly at any boundary of the domain; instead, it is to be determined from the matching 
condition at the singular surface, i.e., the LC (see below). 

Thus, our position differs from that of Beskin (see Beskin et al. 1983, and Beskin 1997) who 
proposed prescribing I (iff) explicitly at the surface of the pulsar. Indeed, if one considers the 
region inside the LC only, then, upon prescribing both iff and I (iff) on the inner boundary 
(the pulsar surface, etc.), one should be able to obtain a solution that is regular at the LC 
by using the regularity condition (20) at the other boundary of this inner region, i.e., at 
the LC. Thus, the solution in this inner region is then completely determined; as such, it is 
totally independent of what happens outside the LC. This does not seem physical; indeed, 
one should then be able to take this solution and continue it smoothly across the LC, thus 
prescribing both iff(l,z) and ^(l,^) at the outer side of the LC. But then the problem of 
finding a solution in the region outside the LC becomes over-determined, as we now have 
two conditions at x — 1, in addition to any conditions at the outer boundary. For example, 
in the very important particular case of the domain under consideration extending all the 
way to radial infinity, one cannot actually prescribe any specific boundary conditions at 
infinity because equation (15) has a regular singularity there. Instead, however, one imposes 
a regularity condition at infinity. In spherical polar coordinates 3 (r, 9, 0) this regularity 
condition can be written as 

H'(iff) 

iff d6 + iff d cot9 = r^oo. (21) 

This regularity condition has a very simple physical meaning: it is the condition of the force 
balance in the ^-direction between the toroidal magnetic field and the poloidal electric 
field Eg (which both become much larger than the poloidal magnetic field B po \ ~ B r as one 
approaches infinity). Thus, in this case the region outside the LC does not have any boundary 
conditions at all (apart from the condition iff = const at the equator), but has two regularity 
conditions, one at the LC, and the other at infinity. These two regularity conditions are 



3 Note that in the rest of the paper we shall use a totally different set of polar coordinates defined in the 
vicinity of the separatrix intersection point, for which we shall use the same notation (r, 9). We hope that 
this does not cause any confusion. 
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already sufficient to uniquely determine the solution ^(x > 1, z) for a given Therefore, 
if, in addition to these two regularity conditions, one also tries to impose the function z) 
along the LC, the system becomes over-determined. Thus we come to the conclusion that 
prescribing I{^>) at the inner boundary is not appropriate. One could, however, try to use 
the outer boundary condition (or the regularity condition at infinity) as the condition that 
fixes To do this, one first sets the inner boundary conditions, picks an initial guess 

for and uses (20) to obtain a regular solution inside the LC; one then continues this 

solution smoothly across the LC and solves equation (15) in the outer region as an initial- 
value problem (with both \I/ and d x ^> specified at x = 1); as a result, one gets a mismatch at 
the outer boundary between the obtained solution and the desired outer boundary condition 
[or, in the case of an infinite domain, one would presumably fail to get a convergent solution 
satisfying the regularity condition (21)]. Then one iterates with respect to until this 

outer-boundary mismatch is zero (or until a solution regular at infinity is achieved). In 
reality, however, such an approach may not be practical as it involves solving an initial-value 
problem for an elliptic equation, which is not a well-posed problem. 

Instead, we advocate the approach adopted by CKF. In this approach one considers the 
regions inside the LC (x < 1) and outside the LC (x > 1) separately. First, one makes an 
initial guess for and prescribes the corresponding boundary conditions for ^ at the 

inner boundary of the the region inside the LC and at the outer boundary of the the region 
outside the LC [or, if this outer boundary is at the radial infinity, imposes the regularity 
condition (21) there]. Then, one uses, again separately in each region, the regularity con- 
dition (20) in lieu of a boundary condition at the LC [in practice, one can think of (20) as 
of a mixed-type Dirichlet-von Neumann boundary condition relating the value of \1/ at the 
surface x — 1 to the value of its derivative normal to this surface]. Thus one obtains two so- 
lutions, one inside and the other outside the LC, that correspond to the same function 
and are both regular at x — 1. However, these solutions also depend on their respective 
boundary or regularity conditions set at the inner/outer boundaries of the domain. Hence, 
in general, although they are both regular at the LC, these solutions are not going to coincide 
at x — 1. The mismatch A^ LC (z) depends on both the chosen function and on the 
inner and outer boundary conditions. One then iterates with respect to in order to 

minimize this mismatch. We believe that for a given set of boundary conditions there should 
be only one choice of for which the mismatch A^lc vanishes and hence the function \1/ 
becomes continuous along the entire LC, i.e., ^i n (x — > l,z) — ty out (x — > l,z). Once this 
special function is found, it then follows from (20), which is satisfied separately on 

both sides of the LC, that d x ty is also continuous, and so the entire solution passes smoothly 
across the LC. Thus, one can say that the function is determined by the regularity 

condition (20) applied separately on both sides of the LC, plus the condition of matching 
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of the two solutions. Another, equivalent way to put it, is to say that is determined 

by the (non-trivial!) requirement that the derivative d x ^(x = 1, z) actually exists at the LC 
[whereby I(^f) is expressed in terms of this derivative via equation (20)]. 

This approach has been implemented successfully in the pioneering numerical work by CKF 
and then subsequently by Ogura & Kojima (2003). Both these groups have used a relaxation 
procedure to arrive at I(^f) that corresponded to a unique solution that was both continuous 
and smooth at the LC. 



3. Behavior near the separatrix Y-point: general considerations 

The main focus of this paper is the behavior of the magnetic field in the vicinity of the point 
of intersection of the separatrix ^ = ^ s between the region of closed field lines (region I) and 
the region of open field lines (region II) with the equator z = (see Fig. 1). We shall generally 
call this point the separatrix intersection point. Our interest in this non-trivial problem is 
fueled by the belief that understanding the key features of this behavior is absolutely crucial 
to devising the proper boundary conditions and providing verification benchmarks for any 
future global numerical investigations of a force-free pulsar magnetosphere. At the same time, 
the treatment of the magnetic field around this very special point is expected to require a 
certain degree of subtlety and delicacy, as noted, for example, by Beskin et al. (1983) and 
by Ogura & Kojima (2003). 

Among the questions that need to be answered are: what is the radial position xq of the 
intersection point (namely, whether this point lies at the LC, Xq = 1, of inside the LC, 
Xq < 1)? and what is the angle 6 at which the separatrix approaches the equator at this 
point, the three a priori possibilities being 8q = (in which case we'll call this point the 
cusp-point), < #o < 7r /2 (the Y-point), and 9 = n/2 (the T-point)? 

In our analysis, we shall assume that equations (1) and (4) that describe our system are 
valid almost everywhere in the vicinity of the separatrix intersection point, i.e., everywhere 
with perhaps the exception of measure-zero regions. Thus, we shall allow for the presence 
in our system of current sheets of infinitesimal thickness, across which the magnetic field 
can experience a finite jump. We limit our consideration to the situation where such current 
sheets can be present only along the separatrix between regions I and II for x < Xq and along 
the equatorial (z = 0) separatrix between the upper and lower open-field regions for x > x . 
Thus, we have one current sheet that lies on the equator at x > x and at x = x splits into 
two symmetrical current sheets (one in each hemisphere) lying along the separatrix ^ = 
(see, e.g., Okamoto 1974). At the same time we assume that our equations apply everywhere 
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else at least in some region around the intersection point, including the portion of this region 
that lies outside the LC (in the case x = 1). 

Let us first assume that the separatrix between closed and open field-line regions reaches 
the LC (i.e., x = 1), and let us consider the separatrix equilibrium condition (e.g., Okamoto 
1974; Lyubarskii 1990). Indeed, whether or not the separatrix contains a current sheet 
(which we assume to be infinitesimally thin), it must satisfy the condition of force balance, 
which is obtained by integrating equation (1) across the separatrix: 4 

(B 2 — E 2 ) 7 = (B 2 — E 2 ) J/ . (22) 

By using equation (12) and also the fact that = in region I, we can rewrite this condition 

as 

[(KM 2 - (KM) 2 ] (1 - - 2 ) = (B 1 /) 2 = ^ , (23) 

n LC x 

where / marks the distance from the intersection point (z = 0) along the separatrix. 
Now, as the separatrix approaches the LC (x — > 1, I — > 0), we get 

2(1 - x) [(Bi ol (l)) 2 ~ (B^(/)) 2 ] = , (24) 

n LC 

where we have defined 

Is=\im (25) 

as the value of the poloidal current of the last open field line (above the separatrix current 
sheet). Since I s is constant along this last open field line, and hence is independent of I, 
we immediately see that in order to satisfy equation (24) in the limit x — > 1, we must have 
either: 

(i) I s = 0, i.e., all poloidal current that flows out of the pulsar must return back to the pulsar 
along open magnetic field lines, with no finite poloidal current flowing along the separatrix; 5 
simultaneously, from equation (24) it then follows that B^ ol (l) = B^(l), and hence no finite 



4 Strictly speaking, equation (1) may not be valid inside a current sheet as other forces may also be 
important in the pressure balance there. However, gravity and the plasma inertia are small because the 
separatrix current layer is assumed to be very thin; and the plasma pressure gradient, while not necessarily 
small inside the current layer, gives, when integrated across the layer, just the difference between the values 
of the pressure on both sides outside the separatrix, both of which are small by assumption. 

5 This scenario is in agreement with the findings by Okamoto 1974 and also with the particular case /3q = 
of the compatibility condition by Bcskin ct al. (1983) and Beskin & Malyshkin (1998). 
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toroidal current flows along the separatrix either (we make a very natural assumption that 
Bp 0l is in the same direction as B^); thus the separatrix in this case is actually not a current 
sheet; 
or, 

(ii) a finite I s 7^ but a divergent poloidal magnetic field in the closed- field region, B^ ol (l) ~ 
1/a/1 — x(l) — > 00, as I — > 0, — > 1. Note that this situation cannot be dismissed 
automatically; indeed, close to the LC the magnitude of the electric field becomes very close 
to that of the poloidal magnetic field, with the infinitely strong electric forces balancing out 
the infinitely strong magnetic ones to a high degree. Then, the separatrix current sheet would 
have to carry back to the star a finite poloidal current I s (thus closing the poloidal current 
circuit), and also a toroidal corotation current with the surface current density /</,(/) divergent 
near the LC: I<j>(l) ~ \Bl \ — ~ — x{l) — > 00; there would also have to be a 

divergent surface charge density on the separatrix: a(l) ~ \E^ ol — ~ 1/ a/1 — x(l) — > 00. 

Note that these conclusions, derived from the separatrix equilibrium condition near the LC, 
are valid regardless of the angle 9q between the separatrix and the equator at the intersection 
point. 

In our analysis in the next section we shall assume that there are no infinitely strong fields 
anywhere in the system, and thus shall dismiss situation (ii) (x = 1, Q = const, finite 
I s 7^ 0) as unphysical. We shall thus concentrate our attention on situation (i) (I s = 0). 
We shall try to analyze the structure of the pulsar equation (15) and obtain the asymptotic 
solution of this equation in the vicinity of the intersection point (x — 1, z — 0). 

4. Behavior near the separatrix Y-point: asymptotic analysis 

In this section we perform an asymptotic analysis of the pulsar equation (15) in the vicinity 
of the separatrix intersection point under the assumption that it is located at the LC, i.e., 
x — 1, Zq — 0. Our approach is actually very similar to the analysis of the magnetic field 
near the endpoint of a non-relativistic reconnecting current layer, performed previously by 
Uzdensky & Kulsrud (1997). 

In addition to assuming x — 1, we shall, in this section, make use of the following two 
assumptions: 

1) the intersection point is a finite-angle Y-point, i.e., the separatrix approaches the equator 
at an angle 9 , < 9 < ir/2, where we measure angles from the radial vector lying on the 
equator and directed toward the star, as shown in Figure 1; 

2) the poloidal current is fully closed in the open-field region, i.e., I s = = ty s ) = 0; this 
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assumption is motivated by the arguments presented in the previous section. 

Furthermore, we also assume symmetry with respect to the equator and therefore consider 
only the upper half-space, z > 0. 

When considering the ^-dependent coefficients in equation (15) in the vicinity of the Y-point, 
we can, to lowest order in \x — 1|, replace x by 1 everywhere except where it appears in a 
combination like 1 — x; thus we can rewrite this equation as 



2{ 



d 2 ^ d 2 ^i „ dty 



+ 



2 UF = - /// W> (26) 



d^ 2 dz 2 \ 0f 

where we introduced a new coordinate, £ = 1 — x. Instead of using coordinates (£,z), it is 
actually more convenient to use polar coordinates (r, 6) defined by 

f = 1- x = rcos9 , (27) 
z = rsind. (28) 

When making the transition to these coordinates, we shall make use of the following expres- 
sions: 



dr\ £ „ fd6\ z sin# . . 

w) = i = ^ e - dr) =4 = — ; (30) 

az/i r \ozJ £ r z r 

as well as the usual expression for two-dimensional (2-D) Laplacian in polar coordinates: 

d 2 m d 2 ^ . 19/ 1 9 2 $ 

+ ^ = A 2 ^ = -- r— +-— -. (31) 



Q£2 q z 2 - r dr\ dr J r 2 d9 2 

In these coordinates the poloidal magnetic field components are written as follows (we do 
not need these expressions now but will need them later): 

1 

B r = p5"nL (32) 
R Z LC rod 

1 

B e = --&-7T- ( 33 ) 
Rlc dr 

Finally, by using expressions (29)-(31), we can rewrite the pulsar equation (26) in our polar 
coordinates (27)-(28) as 

/d 2 ^ 1 d 2 ^\ d^f 1 

(w + *w) ra » 9 + 2 W a * e -7i»* ie = -2 ' (34) 
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Now we need to solve this equation separately in region I (region of closed field lines, < 
9 < 9 ) and in region II (region of open field lines, 9 < 9 < n) and match the two solutions 
together at the separatrix 9 = 9 . Since in this paper we are dealing only with the local 
structure of the magnetic field near the Y-point, we can, without any loss of generality, adopt 
the convention of counting the poloidal magnetic flux \l/ from the separatrix 9 = 9q (instead 
of the usual convention of counting \I> from the rotation axis). Thus, we shall set ty s — 
and, correspondingly, ^ > in region I and \1> < in region II. 

Very close to the Y-point, i.e. at distances much smaller than the Light Cylinder radius (r <C 
1), the system lacks any natural length scale. Hence, we can expect the radial dependence 
of \I/ to be a power law, which enables us to make the following self-similar ansatz: 

tt 7 (r,0) = ^(r,9<9 )=r^f(9), (35) 
*u{r,6) = V(r,9>9 ) = -r a *g(9), (36) 

[we put a minus sign in equation (36) in order to have g{9) > 0.] 

Then, the magnetic field components are given by: 
in region I: 

Bl = * r *- 1 /'^), (37) 
Bl = --Laxr *- 1 /^); (38) 



and in region II: 



r lc 



B? = --i-r^-VW, (39) 
BY = -j^a 2 r^g(9); (40) 



Notice that the condition that magnetic field does not diverge near r = imposes the 
restriction 

«i, «2 > 1, (41) 

whereas for the case where the poloidal magnetic flux function in region I diverges as 1/ y/r 
[i.e., case (it) in § 3], we would have a.\ = 1/2. 

As for the poloidal current function that appears on the RHS of equation (34), it takes 
very different forms in regions I and II. In region I there should be no toroidal field, so 



/(* > 0) = 0. 



(42) 
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In the open-field region II, is not zero but does approach zero in the limit \1/ — > 0. Since 
there is no natural a priori magnetic field scale in the vicinity of the Y-point, we shall again 
employ a self-similar ansatz, i.e., assume that I scales as a power of |^| near ^ = 0: 

/(^ < 0) =g(-^) /3 , (3>0, (43) 

so that the right-hand side of equation (34) becomes 

< 0) = K(-m) 2 P-\ (44) 

where 

K=-/3g 2 >0. (45) 

Upon substituting relationships (35) and (42) (for the closed-field region) and (36) and (44) 
(for the open-field region) into the main equation (34), we obtain two Ordinary Differential 
Equations (ODEs) for the functions f(9) and g(Q), along with a relationship between the 
power-law indices a 2 and (3. In particular, in region I we get a homogeneous linear 2nd-order 
ODE for f(9): 

f"(G) - f'(6) tantf + a^a, + l)/(0) = , (46) 
which is to be supplemented by two boundary conditions: 

f'(6 = 0)=0, /(flo)=0. (47) 

Similarly, in region II we get a non-linear 2nd-order ODE for g(9): 

g"(6) cos 6 - g'(6) sin 6 + a 2 (a 2 + 1) g(0) cos 6 = -k g 1 '^ (6) , (48) 
together with the relationship 

' 8=1 -2^ (49 > 
that follows from the requirement that the LHS of equation (34) scale as the same power 
of r as the RHS, i.e., a 2 — 1 = {2(3 — l)a 2 - [At this point, however, we have to remark that 
one cannot a priori exclude the possibility j3 > 1 — l/2a2; in such a case the poloidal current 
would go to zero near the separatrix so rapidly that its contribution to the pulsar equation 
would become completely negligible. However, as we shall discuss later in this section, it is 
possible to show that no continuous solutions exist in this case.] 

The boundary conditions for equation (48) are 



g(0 o ) =0 = g(9 = n) 



(50) 



- 15 - 



Notice that the constant k that appears on the RHS of equation (48) is in fact unimportant; 
it just sets the scale of variation of the function g{9) and, hence, the overall scale of the 
magnetic field strength. Thus, we can rescale k away by incorporating it into the solution; 
we do this by defining a new variable 

G(0) = K- a *g(e). (51) 

Then, equation (48) becomes 

G"{9) cos0- G'{9) sin# + a 2 {a 2 + 1) G{9) cos9 = - G 1 ~^(9) , (52) 
with the same boundary conditions 

G{6 Q ) = = G(9 = tt) . (53) 

Thus, we see that in our problem the magnetic field structure near the Y-point is completely 
characterized by four finite dimensionless parameters: the separatrix angle 9 and the three 
power-law indices ai, a 2 , and [3. Our goal is to obtain a unique solution of our problem, that 
is to determine the values of these parameters and, simultaneously, determine the functions 
f{9) and g(0). In fact, equation (49) already gives us one relationship between the power 
exponents, but we still need three more relationships to fix all four parameters. Therefore, we 
shall now discuss the conditions that will help us obtain these three additional relationships. 

The second condition [the first being equation (49)] that links our dimensionless parameters 
is rather obvious; it is the condition of force-balance across the separatrix 9 = 9 and it will 
give us a relationship between ol\ and a 2 . We have already discussed this condition in § 3 
[see eqs. (22)-(24)]; with |E| = x|B po i|, > 0) = 0, and I s = 0, this condition can be 
written simply as 

[BUM)Y = [BUM)]" ■ (54) 

Since in our (r, ^-coordinates the separatrix is a radial line 9 = 9 , and since we expect 
the poloidal magnetic field on the two sides of the separatrix to be in the same direction, 
condition (54) simply means that Bj.(r,9 ) = B^(r,9 ). Then, using equations (37)-(40), 
we immediately obtain a second relationship between the power exponents: 

a>2 = on = ot , (55) 

and also a relationship between the normalizations of the functions f(9) and g(9), cast in 
terms of their derivatives at 9 = 9q\ 



f'%) = -9'(0 ) 



(56) 
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The remaining two relationships between the dimensionless parameters come from the prop- 
erties of equations (46) and (52) themselves and cannot, unfortunately, be written out as 
explicit algebraic equations containing these parameters. This fact, however, can prevent us 
neither from describing and discussing the physical and mathematical conditions on which 
these two additional relationships are based, nor from using these conditions to obtain the 
actual unique values of the parameters. 

First, we shall show that there is a unique one-to-one relationship between the separatrix 
angle 9q and the power exponent a. This relationship comes from analyzing equation (46) for 
the closed-field region I. For any given 6 this homogeneous linear equation with boundary 
conditions (47) is an eigen- value problem for the coefficient a = a(a + 1) > 2. It is in fact 
very easy to see why this must be the case. Indeed, the boundary conditions (47) do not give 
us any scale for /; if some function f(9) is a solution of the problem (46)-(47), then Cf(9) 
with an arbitrary multiplier C will also be a valid solution. Thus, the normalization of f{9) 
is arbitrary: if, for given 9 and a, a non-trivial solution of the problem (46) -(47) exists, 
then there will be infinitely many such solutions; in particular, there will be a solution with 
/(0) = 1. Thus, we can impose an additional boundary condition /(0) = 1. Then, however, 
we have three boundary conditions for a 2nd-order differential equation, which means that 
for arbitrary 9q and a the system is over-determined. Hence, such a solution will exist not for 
all values of 9 and a; the condition that it actually does exist gives us a certain relationship 
between 9 and a: for a given 9 we would get an infinite discrete spectrum of the allowed 
values of a. Of these, however, we are interested only in the lowest one because we want a 
solution without direction reversals of the magnetic field, i.e., with / > everywhere. 

From the practical point of view, the easiest way to determine the relationship between a 
and 6*o is the following. We scan over the values of a; for each given a we use our freedom of 
normalization of f{9) to set /(0) = 1. Thus we now have an initial-value problem with two 
conditions at 9 = 0: /(0) = 1 and f'(0) = 0. With these two conditions and with the value 
of a = a(a + 1) given, we can integrate equation (46) forward in 9 (we do it numerically 
using a 2nd-order Runge-Kutta scheme) until the point where f(9) becomes zero. This 
point is then declared the correct value of the separatrix angle 9 for the given a. The 
function 9q{q) we have obtained as a result of this procedure is shown in Figure 2, whereas 
Figure 3 shows the behavior of the function f(9) for several selected values of a. We see 
that 9q(o) is a monotonically decreasing function with 9 (a = 1) = 0.986.. rad [instead of 
6*o(l) = vr/2 that one would get if the Y-point were at a finite distance inside the LCI], and 
9 (a = 2) = 0.614.. rad m 35. 2°. 6 The asymptotic behavior of 9 (a) in the limit a — > oo is 



6 Note that this value agrees with Michel's (1973a) separatrix angle of w 35°, corresponding to the case 
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9o - ¥ n /2a, which is in fact very easy to obtain analytically. 7 Actually, the exact analytical 
fit for our numerical curve 9 (a) appears to be 

^ ) = I^r4rm' (57) 

^ ^a(a + 1) + b 

with b ~ 0.54.... 

With the dependence 9 (a) thus determined, and with j3 and ct2 related to a via equa- 
tions (49) and (55), respectively, we now have everything expressed in terms of a. 

The fourth (and final!) condition that will help us determine a and hence 9 and the rest of 
the power exponents is the Light Cylinder regularity condition. This condition states that 
at the LC 9 = n/2 [which is a regular singular point for equation (52)] the function G(9) 
should be regular, namely, it should be a continuously differentiable function, with finite 1st 
and 2nd derivatives. This condition can be written as 

= c-i(:), m 

and it is just a particular manifestation of the general LC regularity condition (20). We use 
this condition to fix the unique value of a (see the general discussion of this issue at the end 
of § 2). Here is how we do it in practice. 

First, we divide up region II into two sub-regions: region II' (inside the LC: 9 < 9 < n/2) 
and region II" (outside the LC: n/2 < 9 < n). Then, we scan over a; for each value of a, we 
first determine the corresponding value of #0 using the procedure outlined above. Once 9q 
for a given a is found, we solve equation (52) separately in regions II' and II" (again, using 
a numerical shooting method in conjunction with the 2nd-order Runge-Kutta integration 
scheme). In each of these regions we use regularity condition (58) at 9 = n/2 and one of 
the boundary conditions (53) at the other end of the region, i.e., G(9 ) = in region II' and 
G(n) = in region II". Note that the use of the regularity condition (58) guarantees that the 
solutions in each of the two regions are regular. However, in general (i.e., for an arbitrarily 
chosen a), the two solutions obtained in this manner do not match each other at the LC, 
the mismatch AG LC = G n " (n/2) — G H '(n/2) [and hence the mismatch in G'(n/2) related 



of zero poloidal current (see also Beskin et al. 1983); this is because in his analysis he had in fact assumed 
that a = 2 and then derived the corresponding value of 9q. 

7 Indeed, as a — > 00, we expect #0 <C 1, and then we can neglect the /'tan# term in equation (46). As 
a result, we get a simple harmonic equation /"(#) + af(6) = 0, with /'(0) = 0, f(6o) = 0; the solution 
of this equation, positive everywhere in the domain < 9 < # , is f(0) = cos(7T0/2# o ); it corresponds to 
a ~ yfa = n/26 Q . 
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to AGlc via eq. (58)] being a-dependent. We find that there is only one, special value of a 
for which AGlc vanishes and the solution continues smoothly across the LC as one passes 
from region II' into region II". This special value, which we call a , is declared the correct 
value of a. Numerically, we find 

a = a « 1.5045... , (59) 

and, correspondingly, 

P « 0.668... , # ^0.76rad. (60) 

This is the way in which the unique correct solution of the problem is obtained. The functions 
f{9) and G{9) = K~ a g{9) corresponding to this value of a are plotted in Figures 4 and 5, 
respectively. 

Now let us make a little digression and step back to discuss the possibility of finding a 
solution in the case /3 > 1 — XjloLi. Recall that in this case the dominant balance of the 
lowest-order (in r) terms in the pulsar equation does not include the contribution —II'(^>)/2 
due to the poloidal current. At the same time, note that relationships (55) and (56), derived 
from the separatrix force-balance condition, still hold. Hence, the solution in region II would 
have to satisfy a linear equation that is identical to equation (46) for f(9) on the other side 
of the separatrix. This fact, combined with the condition (56), implies that the solution 
just has to pass smoothly across 9 = 9 , and thus this point [which is just an ordinary 
point for equation (46)] does not have any special significance; it is just the point where the 
solution changes sign. Thus, the general solution in this case takes the form \1/ = r a f(9), 
where f(6) satisfies the homogeneous linear equation (46) in the entire domain [0, ir] with 
the homogeneous boundary conditions f'(0) = and /(7r) = 0. In addition, the solution 
has to pass smoothly across the regular singular point 9 = n/2 (the Light Cylinder); the 
regularity condition at this point is simply f'(n /2) = 0. Such a solution would start positive 
with a zero derivative at 9 — 0, then decrease and change sign at some angle 6q < n/2, 
reach a minimum at the LC [since f'(n/2) = 0], and finally increase and go back to zero 
at the equator theta = ir. However, the extra regularity condition at 9 = n/2 makes the 
system over-determined. Indeed, counting two boundary conditions at 9 = 0, ir, an arbitrary 
normalization of the solution [e.g., /(0) = 1], and the LC regularity condition, we now have 
four conditions for our second-order equation with one free parameter a. Thus, unless we 
are somehow extremely lucky, the solution of this problem does not exist. One can show 
that this is indeed the case; it turns out that a solution of equation (46) in the region inside 
the LC ([0, 7r/2]) satisfying both the boundary condition /'(0) = and the LC regularity 
condition f'(n/2) = exists only for even values of a = 2,4,6,..., whereas a solution in 
region ([ir/2,ir]) satisfying both the boundary condition f(ix) = and the LC regularity 
condition exists only for odd values of a — 1,3,5,.... This demonstrates that, in the case 
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P > 1 — l/2a 2 , it is impossible to obtain a solution that is continuous at the LC. This leaves 
us with the case (3=1 — l/2a 2 (for which we have just found a unique suitable solution) as 
the only possibility. 

Thus we have managed to obtain a unique asymptotic solution of the ideal-MHD, force-free 
system (l)-(4) in the vicinity of the separatrix intersection point (x — 1, z — 0). This 
solution, characterized by dimensionless parameters (59) and (60), satisfies the separatrix 
equilibrium condition and is continuous and smooth at the LC. Let us now discuss the 
physical relevance of our solution to the pulsar problem. 

Note that the obtained unique solution of our problem has a finite value of the derivative 
G'(9) at the equator 6 = ir, namely, G'(n) = —2.02, and thus has a non-zero radial magnetic 
field just above (and just below) the equatorial current sheet. This actually spells bad news 
for the applicability of our solution. Indeed, one more condition that needs to be satisfied for 
the solution to be physically relevant and that we have not discussed so far is the condition 

|B| 2 >|E| 2 . (61) 

With |E| related to |B po i| via equation (12), this condition can be expressed as 

E 2 >|B pol | 2 (o; 2 -l). (62) 

We see that this condition is satisfied trivially everywhere inside the LC, but is not auto- 
matically satisfied outside the LC. In particular, in our solution, and in fact in any solution 
characterized by I s = 0, the toroidal field on the last open field line \I/ = is zero, whereas 
the poloidal magnetic (and hence electric) field is not zero, as our explicit solution shows. 
This means that the inequality (62) is violated immediately outside the LC as one moves 
along the equator. More generally, this inequality is violated beyond the so-called Light 
Surface (defined as the surface where |B| = |E|) that emanates from the Y-point. It is in 
fact not difficult to find the location of this Light Surface in our solution. 

Indeed, according to equations (39)-(40), we have 

B 2 pol = ^r*»->[aV(0) + (9'(0)) 2 }, (63) 

n LC 

and, according to equations (19), (43), and (49), 

Bl = 4 (~*) 2/? - ^" ^■ 2a ~ 1 9 2 - l/a (0) • (64) 
Using the result (63), we can rewrite this as 

+ po1 a*G 2 (6) + (G>(6)) 2 • { } 
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Upon comparing this result with 

E 2 - B 2 poX = (x 2 -l) B 2 poX = -2r cos 9 B 2 pol , (66) 

we conclude that the Light Surface B 2 ,=E 2 — B poX emanates from the Y-point at a finite 
angle 9ls which is determined from the algebraic equation [provided that the solution G{9) 
is known]: 

7^- G /a (9 LS ) = - \a 2 G 2 {9 LS ) + (G'(9 LS )) 2 ] cos9 LS . (67) 
la — 1 

For our solution described by the parameters (59)-(60), we find 

# L5 ^2.07rad, (68) 
which corresponds to the angle n — 9ls ~ 62°. 

Beyond the LS (9 > 9 LS ) the force-free equation (1) is not applicable. The solution we have 
obtained in this section is valid in the region 9 < 9 LS , but, because the process of obtaining 
this solution also involved the region 9 > 9ls [e.g., via the boundary condition g(n) = 0], 
our solution may not be the correct solution of the overall problem. 

In fact, we can make an even more general statement: any relevant to our problem finite- 
field magnetic configuration with the separatrix intersecting the equator at the LC cannot 
be almost everywhere (that is, again, everywhere except at a finite number of infinitesimally 
thin current sheets) ideal and force-free beyond the LC! Indeed, as we saw previously, if one 
insists on having the separatrix intersection point at the LC, then one has to contend either 
with having |-B po i| divergent near x = 1 or with having I s = I(^f s ) — 0. The latter case, 
however, is characterized by zero toroidal magnetic field everywhere along the last open field 
line \& = \P S . Provided that the poloidal (radial) magnetic field does not vanish identically 
near the equatorial current sheet, it then follows that E 2 — B 2 oX = (x 2 — l)B 2 oX > B 2 = 0, and 
so the force-free equation (1) becomes inapplicable. And, as our particular solution shows, 
the condition B 2 — E 2 > is violated not only on some singular surfaces of measure zero, 
such as the equatorial current sheet (that would be acceptable to us at this stage, since we do 
expect our simple force- free, ideal-MHD assumptions to break down there anyway), but in a 
volume of non-zero extent in both poloidal directions, e.g., in the region (Q LS < 9 < n,r > 0). 



5. Properties of solutions with the separatrix intersecting the equator inside 

the Light Cylinder 

Thus we have to conclude that if one is to keep the hope of finding a force-free solution that 
would be valid all the way up to infinity (in the limit of vanishing plasma density), or at 
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least to some finite distance beyond the LC, then one has to make concessions regarding 
the anticipated location of the separatrix/equator intersection point. In other words, the 
requirement that B 2 > E 2 outside the LC means that I s must be non-zero and this, in turn, 
forces one to conclude that the intersection point must be located at some finite distance 
inside the LC, i.e., xq < 1. This scenario has actually been the subject of several studies 
in the past 20 years (e.g., Beskin et al. 1983; Lyubarskii 1990; Beskin & Malyshkin 1998). 
According to equation (23), the corresponding finite- I s configuration has to be characterized 
by Bp Ql staying finite in the limit x — > x , and, therefore, the separatrix should approach the 
equator at a right angle, 6 = ir/2, corresponding to a T-point (see Lyubarskii 1990). 

Let us now ask what other characteristic features should such a solution possess? Here 
we present a few simple facts that can be gleaned immediately by inspecting the force- 
free equation (15) and its regularity condition (20). First, since beyond x = xo the last 
open field line, \I> = \I/ S , is assumed to lie along the equator z = 0, the derivative d^/dx 
should be zero everywhere along this line (for x > x ). In particular, we should have 
d^/dx(x = 1, z = 0) = 0, but then it follows from equation (20) that either I(V S ) or I'{^ s ) 
has to be zero. Since here we are considering the case I s ^ 0, we conclude that has 
to approach \J/ S with a zero slope, I'{^ s ) = 0. Thus, we see that ^ cannot be a 

simple linear function, unless it is constant everywhere. Next, by substituting the result 
II'(fy s ) = 0, together with ty(x > xq, z = 0) = = const, into equation (15), we find that 
(1 — x 2 )d 2 ^ /dz 2 (x, z = 0) = for all x > x . This means that the radial magnetic field also 
approaches the equator z = with a zero slope, dB x /dz = 0. 

Now let us see what one can deduce by investigating the regularity condition (20) just above 
the equatorial current sheet. According to what we have just established, above the equator 
the function I(^f) can be expanded as 

J(tt<0)=J a + g (69) 

where f3 > 1 in order to satisfy the condition I'{^>) — * 0, as \I/ — * (here we again set ^ s — 
to simplify notation). Then we find, to lowest order in |^|, 

< 0) = - (3ql s (-mf- 1 . (70) 

Now, the radial magnetic field just above the equator is, in general, not zero, so we can 
expand ^(x, z — > 0) to lowest order in z as 

V(x,z)=zV!(x) + ..., (71) 

where ^i(x) is in fact just the rescaled electric field above the equator, as can be seen from 
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equations (10) and (17): 



xR 2 LC B R (x, z = 0) = -R 2 LC E z (x, z = 0) 



(72) 



Upon substituting equations (70) and (71) into the regularity condition (20), we get 
zV^x = 1) = -lfal a (— \E') /3 ~ 1 = - ^PqlsZ?- 1 [-M^- 1 , 

and, therefore, 

/3 = 2, 

unless = 9 X 9 Z *(1,0) = [and (3 > 2 if = 0]. 



(73) 
(74) 



We can now use this result to analyze the magnetic field in the vicinity of the point (xq, 0). 
Equation (15) in the vicinity of this point can be written as 



;i-*o) 



d 2 ^ d 2 ^ 



dz 2 



x d( 



(75) 



where ( = x — x. Now, just as we did in § 4, we can make a transition to the polar 
coordinates: ( = r cos 9, z = r sin 6. In these polar coordinates equation (75) becomes 



;i-*o) 



rl d 2 ^ 1 d 2 ^ 



+ 
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r <9r <9r 2 r 2 d6 2 . 
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1 + x 2 



X 



— — cost) 
or 



rd9 



sin# 



= -//'(*), (76) 



Again, similarly to what we did in § 4 [see eq. (36)], let us assume that the magnetic flux 
function in region II (the region of open field lines, 9 > 9 = n/2) is a power law of r near 
the point (x , 0): 

#(r,0) = -g(6)r a2 , a 2 > 1 . (77) 



Then, using equations (70) and (74), we can express the RHS of equation (76) as 

RHS = - II' {9) = - 2ql s ^ = 2ql s g(6)r a2 , 
whereas the LHS of this equation can be expressed as 

LHS = -(l-x 2 )r a2 - 2 ia 2 g(6) + a 2 (a 2 -l)g(d)+g"(6)] 

a 2 g{9) cos 6 — g'{6) sin 9 . 



(78) 



1 + X l r a,-l 



x 



(79) 



In this expression, terms of the lowest order in r are those proportional to r° 2 2 ; at the 
same time we see from equation (78) that the RHS of equation (76) scales with r as r a ' 2 . 
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Therefore, we see that in this case of the separatrix intersecting the equator inside the LC, 
x < 1, the contribution of the RHS is completely negligible (as r — > 0). 

Thus, the dominant balance in equation (76) dictates that the terms of order r" 2 " 2 should 
balance each other and hence we get a very simple homogeneous linear ODE for the func- 
tion g(9): 

g"(9) + a 2 2 g(9) = 0. (80) 
The boundary conditions for this equation are 

g(9 = tt/2) = = <7(tt) , (81) 

and so the solution with the lowest a 2 (corresponding to the simplest magnetic field topology 
in this region) obviously is 

g(9) = - g max sin29 , | < 9 < n , (82) 

and, correspondingly, 

a 2 = 2. (83) 

Thus, * n (r, 9) = <7 max r 2 sin2#. 8 

As for region I (the region of closed field lines), the magnetic field there is, to lowest order 
in r, purely vertical and finite at x — > x , so we can write 

ttj(r, 9) = / max C = / max r cos 9 , (84) 
where / max is related to I s via the separatrix equilibrium condition (23): 

/max = 7= • (85) 

a^oV 1 - 

We express hope that these simple findings will be helpful in setting up or checking the 
correctness of future numerical attempts to solve this nontrivial problem. 



8 The fact that we have been able to find such a simple non-divergent solution near the T-point, whereas 
Lyubarskii (1990) has found only a logarithmically-divergent one, can be explained by noting that, as we 
have shown in the beginning of this section, II'(^ S ) must be equal to zero according to the LC regularity 
condition, whereas Lyubarskii assumed that I( V E') ~ W and hence had a finite, non- vanishing II'^s). 
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6. Discussion and conclusions 

In this paper we have considered the axisymmetric force-free magnetosphere of an aligned 
rotating magnetic dipole, under the additional assumptions of ideal MHD and of the unifor- 
mity of the field-line angular velocity [fi(^) = const]. This fundamental model problem is 
of great importance to any attempts to understand the workings of radio-pulsars. 

More specifically we have focussed most of our attention on the structure of the magnetic 
field in the vicinity of the point of intersection of the separatrix (between the closed- and 
open-field regions) and the equator. We call this point the separatrix intersection point. 
The unique singular nature of this point makes it play an extremely important role in the 
overall global problem; in particular, without a thorough understanding of the subtleties of 
the magnetic field behavior near this point, it is impossible to prescribe the correct global 
boundary conditions in any sensible way. 

We start, however, by discussing (in § 2) the basic equations governing the global force- free 
pulsar magnetosphere. We give special attention to the role played by the Light Cylinder 
regularity condition in determining the poloidal electric current. After this general discus- 
sion, the rest of the paper is devoted entirely to the analysis of the separatrix intersection 
point. 

We first consider the separatrix equilibrium condition in the vicinity of this point (see § 3) 
and find that if it lies at the Light Cylinder, then all the poloidal current has to return 
back to the pulsar in the open-field region above the equator, i.e., there should be no singular 
current running in a current sheet along the equator and the separatrix. We then perform 
(in § 4) an asymptotic analysis of the relativistic Grad-Shafranov, or pulsar, equation in the 
vicinity of such an intersection point located at the Light Cylinder. We find a unique self- 
similar solution that can be described by the power law \I/ oc r a , where r is the distance form 
the intersection point and a ~ 1.5045, and by the equator-separatrix angle # 0.76 rad. 
However, a further analysis of this solution in the region outside the Light Cylinder shows 
that a Light Surface (a surface where \B\ = \E\) appears just outside the Light Cylinder; 
in particular, we find that this Light Surface originates right at the intersection point and 
extends outward at a finite angle with respect to the equator. The appearance of a Light 
Surface in this case is, of course, not surprising, taking into account the fact that (in this case 
of the separatrix intersection point lying at the Light Cylinder) the poloidal electric current I 
and hence the toroidal magnetic field have to become zero on the last open field line. We 
therefore conclude that the only possibility for an ideal-MHD force-free magnetosphere above 
the putative equatorial current sheet to extend at least some finite distance beyond the Light 
Cylinder, is for the separatrix intersection point to be located inside, as supposed to right 
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at, the Light Cylinder. (Of course we understand that the exact location of this point can 
only be found as a part of a global solution of the pulsar equation and that it is impossible 
to determine it from our local analysis). 

These findings give us the motivation to consider the case when the intersection point lies 
at some finite distance inside the Light Cylinder. In § 5 we examine the behavior of the 
function J(\&) near the last open field line \I> = ^ s (above the equatorial current sheet) 
for this case; we find that the derivative has to go to zero as \I/ — > \I/ S , whereas the 

current itself can approach a finite value I s 7^ 0. We then perform an asymptotic analysis 
of the magnetic field around the intersection point; it is essentially a somewhat simplified 
and trivialized analogue of our analysis in § 4. We find that the separatrix approaches the 
equator at a right angle, 6 = 7r/2, and that the field in the region of open field lines behaves 
simply as \I>(r, 9 > n/2) ~ — r 2 sin 26*, while in the closed-field region the magnetic field is 
simply vertical and finite near the intersection point (we call this configuration the T-point). 

Finally, we would like to make several remarks regarding the numerical simulations by CKF 
and by Ogura & Kojima (2003): 

1) The midplane boundary conditions (d z ty = 0, R < Rlc\ ^ — const, R > Rlc) adopted 
by both groups have automatically assumed that the separatrix intersection point lies at 
the LC, and thus have precluded them from even considering the possibility xo < 1. It is 
interesting to note, however, that from the magnetic contour plots presented by Ogura & 
Kojima (2003) it does seem that the actual intersection point lies a little bit inside the LC. 

2) Neither CKF, nor Ogura & Kojima (2003) have discussed or even mentioned the separatrix 
equilibrium condition; thus it is not clear whether this condition has been satisfied in their 
simulations. In the light of our present findings and of the fact that both of these groups have 
found I s 0, we suspect that the equilibrium condition has not, in fact, been satisfied in their 
studies, at least close to the separatrix intersection point. This suspicion is strengthened 
by the fact that both groups have reported having experienced some difficulties near the 
separatrix. 

3) Ogura & Kojima (2003) have reported that they had found \E\ > \B\ at some finite 
distance outside the LC, whereas CKF have claimed that in their solution \E\ < \B\ every- 
where. The origin of this discrepancy is not clear. It may be attributed to the difference 
in numerical resolution, although neither of the two groups have conducted very extensive 
convergence studies. 
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Fig. 1. — A schematic drawing of the Y-type separatrix intersection point in the case x — 1. 
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Fig. 2. — The separatrix angle 9 as a function of the power exponent a. 
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Fig. 3. — The function f(9) representing the magnetic flux function in the region of closed- 
field lines (region I) for several selected values of the power exponent a. 




Fig. 4. — The function f{6) representing the magnetic flux function in the region of closed 
field lines (region I) for a = a = 1.50446.... 
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Fig. 5. — The function G(0) representing the magnetic flux function in the region of open 
field lines (region II) for a = ao = 1.50446.... 



